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Significant progress has be e n sade in two areas. A program is 
operational to generate a coordinate system for a two-el snen t airfoil 
with the mesh points concentrated in areas of significant vortfteity, 
l.e. , boundary layer and wake. The development of thenethrd was more 
difficult than anticipated. The "Imbedded" grid net hod developed allows 
a transition from the scale of the main airfoil to the scale of the 
flap. This requirement is essential for the modeling of viscous flows 
over the flap and slat of a multi-element airfoil. Attached to the 
report is a proposal for additional computer tine over the current request 
to explore the fine grid generation problem. 

Progress has also been made in the formulation of the airfoil 
mounted in a 2-D wind tunnel. A Ph.D. thesis has been couplet ed and 
the program is ready for a fine grid and a large number of planes to 
explore the characteristics of a Navler-Stokes solver in a qua si -3D 
case. The program was converted to a fora suitable for the STAR computer. 
Runs will be made to map a 3 dimensional flow field for a wall-airfoil 


intersection with and without lift. 


Mmm I. Mesh Gntntioa for Halt 1-el sweat Airfoil 

The preaemt work duals with the |«Beratioe of bounda ry-f It tod 
coordinate eyttat for • oiltl-eleo^t airfoil by nttodlat the procedure 
developed by Thompson, et al. (It!:. 1). Figures 1 and 2 show a physical 
plane and the corresponding transformed plane In case of two-el caent 
airfoil. H ow e v e r for the convenience of extending the net hod to an air- 
foil having aore than two elements, an alternative arrange men t of trans- 
formed plane that corresp'Mkds to exactly the sane coordinate systen In 
the physical plane Is usad. These ere shown in Figures 3 and 4. 

Now the aaln objective of the further work la to obtain the control 
functions, i.e., the nonhoaogeneous tens of the Poisson's equations, so 
as to get the desired distribution of grid lines In the field. At present, 
two nethods are being investigated for obtaining these control functions. 

The first schane is based on Sorenson's nethod (Ref. 2 and 3). In 
this nethod, the entire region is divided into three sones (in the case 
of two-elenent airfoil} . Then the Grape code is used to generate grids 
for each of the sones. In the Grape code, the specification of spacing 
of the first grid line with respect to boundary, end of angles et which 
£ lines intersect the boundary, are used to obtain the control functions 
p,q,r and a at the top and bottom r\ lines. Thereafter these functions are 
interpolated to obtain the values in the field 

P(C.n) - p(t) e'*^- 15 + r(0 e‘ c(r W “ n) 

Q(C,n) - q(£) e - b * n ~ 1 * + s(C) e^^ax " n) 

Here a, b, c and d ere constants. The continuity across the sonal 
boundary is achieved through these control functions. For example, along 
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a zonal boundary sepai at log two particular rones, it is required that 
the lines, which Intersect that boundary do so at the sane angles with 
respect to boundary (orthogonally In the present case) , with the sane 
•pacing In the direction normal to the boundary. Also, the zonal boundary 
is obtained as a part of the solution, so as to maintain the continuity 
of higher derivatives. 

The second scheme is similar to tb: first one as far as continuity 

across the zonal boundary is concerned. However the control functions 

In the field are obtained by interpolating the control functions calculated 

at all the four boundaries of a zone. The point distribution and the spacing 

of first grid line with respect to boundary are specified at all the boundaries. 

Thus on the top and bottom rj lines, r^ is specified and, therefore r^ 

can be calculated. Using the orthogonality condition at the boundary 

and the spacing t(£), x and y can be calculated as 

ti n 

y £ t(c) 

x * 

n 



x^t(C) 


Thereafter the control functions, P and Q, at the tor' and bottom n lines 
are obtained as 


P 


. Il 


la 


Q - 
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Similarly, on the left end right £ lines, r is specified and 

*n 

therefore r^ can be calculated. Again using the orthogonality condition 
and the spacing t(rj), and y^ can be written as 

y n t(n) 

*{ ' 


'V 


*n t<n) 


'In 1 


and the control functions, P and Q, are obtained using 


r r • r 

P . _ sS 

|r | 2 


Q - - 


r • r 

-n -nn 


In the field, both P and Q can be interpolated as 
F(C, n) - ~ fj F(1 , n) + ~ yj F(l, n) 

+ ^ p(r 1) 4. 1) p(r 1) 

(J - 1) (J - 1) J) 

where 1 < C < I and 1 < n < J . 

Thus a computer program for generating the grids about a multi-element 
airfoil has been developed. At present, two methods for obtaining the 
control functions, ao as to get the desired distribution of grid lines 
in the field, are being investigated. 



Phase II. Numerical Solution of the Time-Dependent, 
Incompressible 3-D Navier-Stokes Equations 

Since the last reporting period, the formulation has been completed 
and run on the CYBER machines. Anutosh Mo it r a completed his Ph.D. thesis 
at this time. He made the conversion of the code to the STAR computer 
to allow execution of a fine mesh and multiple planes to make a 
meaningful simulation of the flow at large Reynolds numbers. This 
transition modification is now complete and the long runs will be executed 
soon. The initial runs will be at zero lift. Then the airfoil was placed 
at an angle of attack to develop flows of interest at the airfoil-wall 
interface. Finally an aifoil with a plain flap deflected to a large 
angle will be run to explore separation phenomena. The grid mesh will 
be modified to concentrate the mesh points in the wake regions. 
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Figure 1. Physical Plane 



Figure 2. Transformed Plane. 







